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Abstract 



We report high spectral resolution observations of the CO vibrational overtone emission from 
the young stellar object 1548C27; our observations include both the v=2-0 and v=5-3 band- 
head regions. These data and similar observations of the young stellar object WL16, reported 
in a previous contribution to this journal (Letters) , provide some of the most compelling ev- 
idence to date for the existence of inner disks around young stars. We describe the simple 
procedure that we use to synthesize bandhead emission from disks including the effect of ther- 
mal dissociation of CO and non-LTE excitation of the vibrational levels. Using this spectral 
synthesis procedure to extract the kinematics and physical properties of the emitting gas from 
the overtone data, we show how these high signal-to-noise data are also powerful probes of the 
stellar and inner disk properties of these systems. Our modeling is consistent with the iden- 
tification of WL 16 and 1548C27 as Herbig AeBe stars with stellar masses of approximately 
2 and 4 M©, respectively. Thus, the kinematic signature of rotating disks in the overtone 
spectra of these sources provide strong support for the role of accretion disks in the formation 
of intermediate mass stars. For both WL 16 and 1548C27, we interpret our modeling results as 
indicating that the overtone emission arises from a temperature inversion region in the inner 
disk atmosphere. We also find evidence for suprathermally broadened lines and are able to 
place useful constraints on the radial temperature and column density distributions of the CO 
line-formation region of the disk atmosphere. Given these deduced properties, we discuss the 
constraints that our observations place on the physical processes responsible for the overtone 
emission in these sources. 

Subject headings: accretion, accretion disks — infrared: stars — line: profiles — line: forma- 
tion — stars: circumstellar matter — stars: formation — stars: pre-main sequence 

1. Introduction 

Circumstellar disks play a fundamental role in the formation of stars and planets. Stars are 
believed to build up a substantial fraction of their main-sequence mass by accretion through 
circumstellar disks and these disks constitute the reservoirs of mass for the formation of both 
planets and close stellar companions. Various pieces of observational evidence assembled over 
the past decade have been used to argue for the presence of accretion disks around young 
stars. These include the spectral energy distributions of young stars (e.g., Adams, Lada, & 
Shu 1987), the asymmetry of their optical forbidden line profiles (e.g., Edwards et al. 1987), 
their millimeter and submillimeter line and continuum emission (cf. Ohashi & Hayashi 1995; 
Koerner & Sargent 1995 and references therein; Keene & Masson 1990; Lay et al. 1994), 
and the cross-correlation linewidths of optical and near-infrared lines of FU Ori objects (e.g., 
Hartmann & Kenyon 1987a, b). More recently, direct kinematic evidence for the existence of 
circumstellar disks has been obtained from high spectral resolution observations of the CO 
overtone emission of young stars (Carr et al. 1993, Chandler et al. 1993, 1995). Since the 
CO overtone lines arise from the inner regions of these disks, detailed studies of this emission 
can be used to determine the properties of disks within several tens of stellar radii of young 
stars. 

This inner disk region (< 1 AU) is of particular interest as it is believed to be dynamically 
active, giving rise to energetic inflows and outflows of gas which remove substantial amounts 
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of mass and angular momentum from the system and regulate the rotation periods of young 
stars (Shu et al. 1994; cf. Edwards et al. 1993; and Bouvier et al. 1993). Although regions of 
this size cannot be imaged directly with present observational techniques, they can be studied 
with high resolution near-infrared spectroscopy of lines that selectively probe the physical 
conditions of the region. The exceptional stability of the CO molecule makes it a versatile 
probe of this region and the excitation temperatures and critical densities of the overtone 
bands near 2.35/xm are particularly well suited to the study of the high densities (up to 10 15 
cm -3 ) and temperatures (1000-6000K) that characterize the inner regions of young stellar 
objects (YSOs), e.g., as first demonstrated by Scoville et al. (1983) for the BN object. 

With the advent of a new generation of sensitive near-infrared spectrometers (e.g., CSHELL, 
Tokunaga et al. 1990), it is now possible to study low- and intermediate-mass YSOs in the 
near-infrared at high spectral resolution. Using these spectrometers, even faint, low-luminosity 
(~ IOLq) YSOs can now be observed with a resolution up to 40,000 and at a signal-to- 
noise that was previously only feasible for bright, high-luminosity sources (e.g., Mitchell et al. 
1989, 1990) and FU Ori objects (Hartmann & Kenyon 1987a, b). We have previously reported 
spectroscopic observations of the v=2-0 bandhead region of the embedded sources SVS 13 
(Carr & Tokunaga 1992) and WL 16 (Carr et al. 1993). Similar observations have been 
reported by Chandler et al. (1993) for several additional sources. Our recent observations of 
the v=2-0 and v=5-3 bandhead regions of the Herbig AeBe star 1548C27 are presented here. 
A striking aspect of these studies is that the bandhead often shows the characteristic shape 
of emission from a rotating disk. The WL 16 and 1548C27 spectra provide particularly good 
examples and clearly demonstrate the existence of inner disks in these systems. 

Winds from low-luminosity YSOs may also produce substantial CO overtone emission (Carr 
1989; Chandler et al. 1995; Saner et al. 1995). Calculations of the thermal and chemical 
structure of winds from these sources show that CO is a robust chemical product and that the 
inner wind region is warm enough to excite overtone emission (~ 2000K; Ruden, Glassgold, 
& Shu 1990). Since the densities of winds rapidly drop below the critical density for the 
overtone bands (cf. Appendix B), a significant wind contribution to the overtone emission is 
expected to be restricted to YSOs with large mass loss rates. It is less clear whether a large CO 
abundance can be expected to survive possible photodissociation by the stellar radiation field 
of the earlier stellar spectral type sources in the present study. In any case, overtone emission 
from winds is expected to show a line asymmetry due to the occultation of the receding flow 
by optically thick inner disks. Since the wind is likely to be cooler than the star and to have a 
temperature similar to that of the inner disk region (Ruden, Glassgold, & Shu 1990; Chandler 
et al. 1995), the absorption of stellar and inner disk continuum photons by the wind will 
introduce an additional line asymmetry. The absence of these asymmetries in the overtone 
spectra in the present study suggests that the emission from these sources is not dominated 
by the contribution from a wind. 

In this paper, we demonstrate that our overtone observations of WL16 and 1548C27 are well- 
explained by a disk interpretation for the emission. We develop and implement a simple 
spectral synthesis technique in order to illustrate the constraints that these high resolution, 
high signal-to-noise data place on both basic stellar properties (e.g., masses and radii) as 
well as detailed disk properties (e.g., radial variations of temperature, column density, and 
intrinsic linewidth). In previous publications, we used a spectral synthesis procedure developed 
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by Carr et al. (1993), which generalized previous work addressing isolated lines (e.g., Huang 
1972, Smak 1984) to the bandhead. The synthesis problem was formulated in terms of an 
effective emissivity as a function of disk radius; inner and outer disk radii were introduced as 
fitting parameters. In addition, the line emitting region was assumed to correspond to the 
entire vertical disk column density and to be optically thin in the continuum (see also Carr 
1989, Chandler et al. 1995). 

In the present application, we develop a more general synthesis method formulated in terms of 
the physical properties of the CO line-formation region of the disk atmosphere. We use a two- 
layer approximation to this atmospheric region: the continuum forms in the lower layer and 
the CO lines form in the upper layer. Thus, the two-layer model approximates the situation 
of an optically thick disk with a temperature inversion in the upper disk atmosphere (cf. 
Calvet et al. 1991). We adopt radial power-law distributions of temperature in both layers 
and of column density in the upper layer. Our analysis then yields the radial variation of 
physical conditions that must be explained by a more fundamental future theory. We avoid 
the introduction of arbitrary inner and outer disk radii by explicitly considering the abundance 
and excitation of CO. In addition to placing tight constraints on the temperature and column 
density distributions of the emitting layer, we find that the intrinsic (local) linewidth of the 
emitting gas is significantly broader than thermal Doppler values, probably symptomatic of 
supersonic MHD turbulence. 

We begin, in §2, by defining the two-layer disk atmosphere model and outlining the procedure 
used to synthesize the CO spectrum of a rotating disk from known level populations. In the 
next three sections, we discuss the applicability of LTE level populations and the circumstances 
under which departures from LTE become important. In §3, we model the v=2-0 bandhead 
spectrum of WL 16 assuming LTE in order to illustrate the parameter dependence of the two- 
layer model and extract the basic physical properties of the emitting gas. In §4, we develop 
a non-LTE theory for the level populations and draw attention to the diagnostic power of 
jointly fitting both the v=2-0 and v=5-3 bandhead regions. The utility of modeling these two 
spectral regions is demonstrated in §5 where we present our observations of the Herbig AeBe 
star 1548C27 in these spectral regions and use our spectral synthesis procedure to deduce 
the physical properties of the CO line formation region in this source. Given these deduced 
properties, in §6 we discuss the evolutionary status of WL 16 and 1548C27 and comment on 
the physical processes that are likely to be responsible for the CO emission in these sources. 

2. A Two-Layer Disk Atmosphere Model 

We model the CO overtone emission as arising from the inner ( ^ 1AU) region of a geomet- 
rically thin disk in Keplerian rotation about a star of mass M*. We treat the atmosphere 
as two adjacent layers, each of which is vertically homogeneous. The overtone lines form in 
the upper ("line") layer, which extends vertically over a thermal scale height and has negligi- 
ble continuum opacity. The latter assumption is checked a posteriori by computing the H~ 
continuous opacity. There is no scattering in the upper layer and the only line opacities are 
those of the CO lines themselves. The lower ("continuum") layer is optically thick in both 
the overtone lines and in the 2.3/xm continuum. 

We adopt power-law radial temperature distributions for the continuum and line layers: 
T c {r) = T c o(r /ro)~ p<: and Tt(r) = Tio(r /ro)~ Pl where ro is a fiducial disk radius. The line 
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layer is additionally described by a vertical mass column density distribution that also varies 
as a power-law with radius: (Tt{r) = &to(r/ro )~ qi . We assume solar abundances of gaseous 
carbon and oxygen in calculating the chemical equilibrium abundance of CO in this layer. The 
additional parameters required to synthesize the absolute strength and shape of the spectrum 
near the CO bandheads are: the local line profile function (<f> v ), the distance and A'-band 
extinction to the system (<f, Ak), the inclination of the rotation axis of the disk to the line of 
sight (z), and the stellar luminosity and radius R*). The last two quantities determine the 
stellar contribution to the 2.3 /im continuum which is assumed to be blackbody. 1 Of these, the 
stellar and system parameters — M*, L#, i?*, <f, Ak, and i — are relatively well constrained by 
the gross properties of the CO overtone emission and other observational data (see, e.g., Carr 
et al. 1993). The disk properties — T c , Tt, at, and <f> v — which are the quantities of primary 
interest here, are well constrained by the type of detailed modeling that we discuss below. 

The synthesis procedure allows the introduction of inner and outer disk truncation radii (i?i n , 
-Rout) as model parameters and explicitly considers the variation of CO abundance and exci- 
tation with disk radius. Both disk truncation and these latter processes may limit the radial 
extent of the CO emission. For example, disks may be physically truncated from the inside 
at distances of several stellar radii by stellar magnetospheres, as suggested by a variety of 
recent observational and theoretical results (c.f. Bertout et al. 1988, Konigl 1991, Bouvier 
et al. 1993, Edwards et al. 1993, Hartmann et al. 1994, Shu et al. 1994). In addition, the 
dissociation of CO at temperatures Tt ^ 5000 K will naturally introduce an effective inner 
radius to the CO emission; in contrast to the case of physical disk truncation, contributions 
to the 2.3//m continuum may arise from within the CO dissociation radius. Disks may also 
be physically truncated from the outside by the formation of giant planets or the presence of 
nearby stellar companions (e.g., Lin & Papaloizou 1993, Mathieu 1994). Since the YSOs in 
the present study are not close binary systems and giant planets are believed to form outside 
the water-ice condensation radius (at T ~ 150 K) — well beyond the region responsible for 
the CO overtone emission — we generally set the outer radius of the disk i? ut at infinity. In 
this case, effective outer radii for the CO emission are produced by either dust formation or 
departures from LTE populations (see §§4 and 6). 

To synthesize the bandhead emission, the disk surface is first divided radially and azimuthally 
into finite elements. For each element, we calculate the emergent intensity 

I v = B v (T c )e- T »l» + 5,(1 - e"^), (1) 

and then add a frequency shift appropriate to the projected rotational velocity of the element. 
In the equation above, B V (T C ) is the Planck function at the temperature of the continuum 
layer, // = cos z, and S v and t v are the source function and total vertical optical depth through 
the upper layer of the element. Both the source function and optical depth of the upper layer 
may include contributions from discrete lines that overlap in frequency. The line layer is 
assumed to be spatially thin so that the properties of the disk atmosphere are approximately 
constant along the inclined ray path. Having calculated the intensity from each element, the 

We ignore the stellar contribution to the CO overtone lines. While this approximation is 
suitable for the probable early-type sources discussed in this paper, the contribution of the 
stellar absorption lines to the composite YSO spectrum is required to model later-type sources 
such as DG Tau (cf. Chandler et al. 1993). 
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contributions of all elements not shadowed by the star are then summed to obtain the disk 
contribution to the spectrum. Finally, we add the stellar (blackbody continuum) contribution 
to the spectrum, accounting for the shadowing of the star by the disk. The resulting spectrum 
is extincted by the adopted value of Ak to produce the final synthesized spectrum. 

3. LTE Model 

In order to illustrate the dependence of the bandhead emission on the parameters of the two- 
layer model, in this section we apply the above spectral synthesis procedure to our u=2-0 
data for WL 16 assuming LTE level populations. The rotational and vibrational constants, 
and the energy levels are taken from Mantz et al. (1975). The use of LTE level populations 
allows us to extract the basic kinematics of the emitting gas from the bandhead spectrum in 
a simple way, i.e., without having to address the additional complication of non-LTE effects 
which generally bear weakly on the deduced kinematics. The LTE results also form a basis 
of comparison against which we can measure the magnitude of the non-LTE corrections that 
are computed in §§4 and 5. 

3.1. Spectral Synthesis of the v=2-0 Bandhead Region 

The spectral region covered by our observations of the u=2-0 bandhead emission from WL 16 
(cf. Carr et al. 1993) includes the R39-R62 lines which approach and retreat from the bandhead 
located at a rest wavelength of 2.29353 /im (Fig. 1). As described by Carr et al. (1993), the 
bandhead emission from WL 16 displays the characteristic shape of emission from a rotating 
disk: a blue wing, a shoulder, and an intensity peak redward of the rest wavelength of the 
bandhead. That this is the general shape expected for bandhead emission from a rotating disk 
can be demonstrated by imagining the convolution of the profile of an isolated line from a 
rotating disk with the distribution of the lines near the u=2-0 bandhead (see also Carr 1995). 

This convolution procedure is shown schematically in Figure 2 where the double-horned ro- 
tational broadening function (heavy line) is the profile of an isolated line from an inclined 
Keplerian disk that has a monotonically decreasing intensity distribution between an inner 
and outer radius. The rest spectrum of lines near the bandhead (light line) includes line blend- 
ing due to local line broadening sources and assumes LTE level populations at a typical inner 
disk temperature. The blue wing is built up as the wing of the isolated line profile moves past 
the bandhead (Fig. 2a). A shoulder is formed when one of the horns of the isolated line profile 
reaches the bandhead (Fig. 2b). As the horn moves past the bandhead it encounters other 
lines longward of the bandhead; the emission consequently remains approximately constant, 
producing the shoulder. The red peak forms when the other horn of the isolated line profile 
reaches the bandhead (Fig. 2c). When both horns have moved past the bandhead, the overlap 
of the isolated line profile with the lines redward of the bandhead produces the slow descent 
redward of the red peak (Fig. 2d). In the observed bandhead spectrum, the velocity separation 
between the blue shoulder and the red peak is approximately twice the projected rotational 
velocity at the outermost radius of the CO emission, and the maximum observed velocity in 
the blue wing is the projected velocity of CO at the innermost radius of the emission. 

The convolution described above is essentially the procedure used by Carr et al. (1993) to 
estimate the stellar mass and radius of WL 16. While suitable for this purpose, this procedure 
is approximate in that the relative strengths of the lines near the bandhead are assumed to 
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be independent of disk radius. The dependence of the emergent spectrum on radial variations 
of temperature and column density enters only in the rotational broadening function which 
is effectively the isolated profile of an "average" line at the bandhead. The present modeling 
procedure removes both of these approximations and allows us to place physical constraints 
on the properties of the inner disk region. 



WL 16 is a low luminosity YSO in the p Oph cloud (<f=160 pc). Although originally classified 
as an embedded (Class I) object (Wilking, Lada, & Young 1989), the measurement of strong 
3.3 /xm aromatic hydrocarbon emission (Tokunaga et al. 1991) and a preliminary analysis of 
our observations of the v=2-0 bandhead emission (Carr et al. 1993) imply an earlier spectral 
type and higher mass for WL 16 than is usually ascribed to Class I objects (e.g., Adams et 
al. 1987). The bolometric luminosity of WL 16 has been estimated as 18 — 22 Lq (Wilking et 
al. 1989; Cabrit & Andre 1991). Carr et al. (1993) estimate the extinction to WL 16 to be 
Ak = 2 based on the observed ( J-H) color of WL 16 and assuming that its intrinsic color is 
similar to that of very active pre-main-sequence stellar systems. 

The analysis by Carr et al. (1993) of the u=2-0 bandhead emission revealed that parameters 
of M* ~ 2.5 M©, i?* < i?j n ~ 5i?0, and i ~ 60° are required to fit both the observed velocity 
extent and flux of the bandhead emission. As discussed by these authors, an inclination of 
i = 60° represents a compromise between the projected emitting area required to fit the 
bandhead flux and the projected rotational velocities required to fit the velocity extent of the 
emission. Several general considerations constrain the stellar mass to be ~ 2.5 Mq given the 
assumed extinction to the source. On the one hand, the stellar mass cannot much exceed 
~ 2.5 Mq for any pre-main-sequence age given the bolometric luminosity of the system. On 
the other hand, a smaller stellar mass would require that the emission arise closer to the star 
to maintain the velocity extent of the bandhead; an origin at smaller radii will tend to reduce 
the strength of the emission since, for a given bandhead shape, the emission scales as R'm . 
The impact of the reduced emitting area on the bandhead flux can be offset only to a limited 
extent by adopting a higher temperature for the emitting gas since collisional dissociation of 
CO becomes important above temperatures of ~ 5000K. 

To illustrate the additional information that can be extracted from the WL16 v=2-0 data with 
the present spectral synthesis method, we reanalyze the data using the two-layer LTE model 
described above. We assume that nearly all of the system luminosity arises from the star and 
adopt d = 160 pc, i = 60°, M* = 2.5 M©, Ak = 2, and = 22 Lq; the remaining parameters 
which describe the properties of the disk — T c o, p C} Tio } pi } cr^o, qi } and <f> v — are varied to fit 
the observed bandhead spectrum. For the continuum layer, we adopt a radial temperature 
dependence in which p c = 3/4, as expected for an active viscous or passive reprocessing disk. 
The adopted disk continuum temperature corresponds to a disk luminosity of 



3.2. Application to WL 16 




in 




and implies a disk mass accretion rate which can be obtained from the relation 




L D -F(1-A)L* 



(3) 
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where the second term accounts for passive reprocessing of incident stellar radiation (Adams 
& Shu 1986). The quantity T is the fraction of the stellar luminosity that is intercepted by 
the disk {T = 1/4 if R- m = i?*) and A is the effective scattering albedo. Values for the other 
disk parameters T c o, T^o, and a to are chosen to reproduce both the strength of the 2.3/xm 
continuum and the contrast of the bandhead above the continuum. More detailed spectral 
features such as the shape of the blue wing, the sharpness of the shoulder, the location of the 
red peak, and the spectral shape redward of the peak constrain pt and qt- 

With these assumptions, we find that we can fit the strength of the 2.3/xm continuum if the 
M* = 2.5 M©, L* = 22 Lq star lies near the stellar birthline (i?* = 5.0 Rq } Palla & Stahler 
1993) and the continuum layer has a temperature of T c o = 3000K at ro = 5Rq. In order to 
fit the velocity extent and strength of the blue wing, Tt must be ^ 5000K at 5Rq. These 
high temperatures correspond to large line source functions without appreciable dissociation 
of CO, and thereby maximize the bandhead flux from these radii. The location of the shoulder 
and red peak indicates that most of emission arises from within 30 Rq; the sharpness of these 
features indicates that the emission must decrease relatively sharply beyond this radius. If 
the rapid decrease in the emission beyond 3Oi?0 is due, in part, to the reduced sensitivity 
of the v=2-0 bandhead lines to gas below 1500K (see below), the constraints on the line 
layer temperature are roughly Tt ~ 1500 K at ~ 30 Rq and Tt ~ 5000 K at ~ 5 Rq which 
corresponds to a temperature profile with Tto ~ 5000 and pt — 3/4. Flatter temperature 
profiles that roughly preserve the locations of the shoulder and red peak and velocity extent 
of the the blue wing have lower temperatures within ~ 3Oi?0 and require somewhat steeper 
column density distributions to fit the sharpness of the shoulder and red peak. 

Given these restrictive requirements, it is very difficult to fit the strength and shape of the 
bandhead with thermal Doppler profiles. Within the context of the two-layer LTE model, 
we can obtain a good fit only within a very restricted range of parameters about oto = 
5000 gem" 2 , qt = 4.5, Tto = 6000K, p t = 3/4, and i? in = i?* (Model L12; see Fig. 3a 
and Table 1). Since the wings of gaussian profiles are weak, it is difficult to emit much 
flux between the line cores without these large column densities and the higher line layer 
temperature. These column densities are, in fact, so large as to violate our assumption of a 
negligible 2.3/xm continuum optical depth r in the line layer: the opacity due to H~ alone 
corresponds to r = 700 at l.li?* and r > 1 within 1.8i?*. 

Including this continuum optical depth in the model would significantly worsen the fit. In 
Figure 3a, we are able to fit the strength and extent of the blue wing and red peak because 
the high column density close to the star yields a sufficient abundance of CO to keep the 
optical depth of the overtone lines appreciable even at Tt = 6000K. If we were to include 
the line layer continuum optical depth in the model, the high optical depth close to the star 
would significantly reduce the contrast of the line emission from this region and result in an 
excessively large 2.3/xm continuum as well as both an inadequate velocity extent and overall 
strength for the bandhead emission (cf. Fig. 3d). In general, somewhat smaller line layer 
column densities can produce the same bandhead flux if the line layer represents the entire 
disk column density within a few stellar radii, i.e., if the continuum layer is absent. In the 
present case, however, even the absence of a continuum layer cannot sufficiently reduce the 
required column density to values corresponding to reasonable H~ opacities within ~ 1.6R*. 
We conclude that, given the adopted extinction to the source, thermally broadened lines 
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cannot account for the CO emission from WL16. 

Alternatively, line layers with column densities of ~ 100 gem -2 at IORq and underlying opti- 
cally thick continuum layers can produce good fits to the strength of both the bandhead and 
the 2.3/xm continuum if the intrinsic line profiles are broader than thermal. With Lorentzian 
lines, these fits require systematically lower column densities than in the case of thermal 
Doppler lines because more flux can be emitted in the stronger wings of Lorentzian lines. 
Figure 3b (Model L14) shows the fit obtained with LTE populations, Lorentzian line profiles 
with a constant line width (half- width at half maximum) of 7 = 0.5 kms -1 , aio = 370 gem -2 , 
qi = 3.8, pi = 0.78, and all other parameters the same as in Model L12 (cf. Table 1). At the 
lower line layer column densities of this model, the CO dissociation radius lies at a somewhat 
larger distance than in Model L12: in Model L14 the abundance and optical depth in the 
overtone lines become substantial beyond 1.2R#. Beyond this distance, the continuum optical 
depth due to H - is negligible, consistent with our assumptions. We are able to fit the extent 
of the blue wing despite the larger dissociation radius because of the stronger wing of the 
Lorentzian line. With these parameters, the shadowing of the star by the disk reduces the 
stellar contribution to the 2.3/xm continuum to ~ 0.29 Jy. The star, in turn, shadows the part 
of the disk that contributes line emission at low projected velocities, enhancing the sharpness 
of the shoulder. The (partially shadowed) disk continuum layer contributes an additional 
~ 0.12 Jy to the 2.3/xm continuum and ~ 2Lq to the system luminosity (eq. 2). 

Because the lines in this model are much less optically thick than in the previous case, there 
is a larger region of parameter space within which we can find acceptable model fits. For 
example, a similar fit is possible with a column density distribution as flat as qi = 2 (with 
aio = 90 gem -2 , Tio = 5500K, and pi = 3/4). Temperature distributions as flat as pt = 0.5 
can also fit the data if the column density distribution is steeper (e.g., aio = 400 gem -2 , 
qt = 4.5, and Tio = 4700K). Larger Lorentzian linewidths further decrease the column density 
required to produce the strength of the bandhead and tend to flatten the slope of the spectrum 
redward of the peak. Using somewhat steeper temperature gradients to compensate for this 
flattening will produce a good fit. For example, for 7 = 2.5 kms -1 , adopting Tio = 5900K, 
pi = 0.78, aio = 75 gem -2 , and qi = 3 will also fit the data. 

WL 16 is unusual among low luminosity YSOs in showing strong aromatic hydrocarbon emis- 
sion features (e.g., Hanner et al. 1992), a property that has been interpreted as evidence for an 
early stellar spectral type. High resolution if-band spectroscopy of WL 16 is also consistent 
with an early spectral type (see §3.4). The system luminosity of 22Lq and deduced stellar 
mass restrict the evolutionary state of such a star. In order to investigate the possibility of 
an earlier stellar spectral type for WL 16, we consider whether WL 16 is an A-star with mass 
M* = 2.5M0, radius i?* = 2.O8i?0, and luminosity = 22Lq. Figure 3c (Model L15) shows 
the synthesized spectrum obtained with these stellar parameters, Lorentzian line profiles, and 
the following set of disk parameters: ro = 5Rq } T c o = 3000K, Tio = 581 OK, pi = 0.78, 
aio = 310 gem -2 , qi = 3.7, and 7 = 0.5 kms -1 (cf. Table 1). The star contributes ~ 0.1 Jy 
of the continuum; the disk extends up to the surface of the star and contributes ~ 0.2 Jy to 
the 2.3/xm continuum and ~ ALq to the system luminosity. The (inner) dissociation radius 
is effectively ~ 5 Rq } as in Model L14. It would be difficult to distinguish between this disso- 
ciation radius and a disk truncation radius i?j n ;$ 5 Rq because the additional (unshadowed) 
stellar contribution to the continuum in the latter case is comparable to the contribution to 
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the continuum from the truncated disk region. Although the fit to the bandhead emission is 
comparable to that obtained with Model L14, the continuum is reduced by a third due to the 
earlier stellar spectral type. 2 

3.3. Discussion of WL 16 Results 

The detailed modeling procedure presented above confirms the result obtained by Carr et al. 
(1993) that the parameters M* = 2.5 M©, i?* < i?j n = 5 Rq } and i = 60° provide a good fit to 
both the observed velocity extent and flux of the bandhead emission. We are able to explain 
the inner disk radius of 5 Rq as either a CO dissociation radius (Model L15) or a physical 
truncation radius. The disk may either extend up to the surface of a 5Rq star on the stellar 
birthline (Model L14) or it may be truncated due to the presence of a stellar magnetosphere 
that extends to a similar distance (a variation of Model L15). 

By fitting the detailed shape of the bandhead emission, e.g., the location and sharpness of the 
shoulder, we find that that the v=2-0 emission must be truncated relatively rapidly beyond 
3Oi?0. This effect can arise for a variety of reasons and probably does not indicate a physical 
outer edge to the disk. For example, at temperatures below 1500K, the lines that make up 
the v=2-0 bandhead are not well excited even if the level populations are in LTE. In the 
models shown, a steep column density gradient is used to emphasize this outer radius to the 
emission. Flatter column density distributions are possible if dust is present in the upper disk 
atmosphere (i.e., the line layer) below 1500K. Because the CO lines and the 2.3//m continuum 
would then form in the same layer this would eliminate the contrast of the CO emission above 
the continuum beyond this radius, creating an effective outer radius to the emission. The role 
of non-LTE effects in creating effective outer radii to the overtone emission are discussed in 
the next section. 

We find that suprathermal line widths are needed to fit the strength and shape of the bandhead 
emission. The bandhead emission is well fit with Lorentzian lines characterized by 7 = 0.5 — 
2.5kms _1 . In §6 we offer a physical interpretation for these enhanced line widths. We also 
find that the required line layer column densities (~ 100 gem -2 at 10 Rq) are much lower 
than the values (~ 10 5 gem -2 ) expected by extrapolating models of the "minimum solar 
nebula" to these radii. We interpret this result as indicating that the line layer represents a 
temperature inversion region in the upper atmosphere of a disk that has a much larger total 
column density, and note that this result also indicates the sensitivity of the CO overtone lines 
to relatively small amounts of circumstellar matter. Similar column densities may be present 
in the continuum "gaps" or "inner holes" of YSO disks (cf. Strom et al. 1989; Hillenbrand et 
al. 1992; Lada & Adams 1992; Marsh & Mahoney 1992; Mathieu 1994). 

The disk continuum temperature used in all fits corresponds to a disk mass accretion rate 
that is consistent with zero but is in all cases ^ 2 X 10 -7 Mq yr -1 . The 2Lq radiated by the 
disk in Model L14 can be entirely attributed to reprocessed starlight if the effective albedo is 
0.68 (eq. 3 with T = 1/4). Similarly, the ALq disk luminosity in Model L15 could be due to 
reprocessed starlight alone if the effective disk albedo is 0.2. Alternatively, the disk luminosity 



2 Note that some fraction of the 2.3//m continuum has been added "by hand" to produce the 
synthesized spectra shown in Figures 3a,b,c (Models L12, L14, L15). The discrepancy is minor 
in the case of Models L12 and L14 but accounts for ~ 1/3 of the continuum in Model L15. 
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in both models can be attributed to energy release by disk accretion at ~ 2 X 10 Mq yr _1 
if the effective albedo is unity. This accretion rate is low compared to typical values for both 
Class I objects (e.g., Adams et al. 1987) and Group II Herbig Ae stars (Hillenbrand et al. 1992), 
consistent with the lack of evidence for either an outflow from WL 16 or other phenomena 
commonly associated with active Class I sources such as SVS 13 and L1551 IRS5. 

The synthetic spectra shown in this section indicate the degree to which modeling the v=2-0 
bandhead spectral region assuming LTE level populations allows us to constrain the physical 
properties of the disk. We find that there is some trade off between the line layer temperature 
and column density profiles required to obtain a good fit. This non-uniqueness is due to both 
the assumption of LTE and the choice of the v=2-0 bandhead as the region to model. Since the 
assumption of LTE fixes the source function, this leaves the optical depth as the only variable 
in the radiative transfer that is affected by a given temperature and column density. Since the 
lines in the v=2-0 bandhead region are excited at similar temperatures and moreover probe 
a restricted range in optical depth at a given temperature and column density, the trade off 
between temperature and column density in the optical depth will work almost equally well 
for all the lines in the spectral region, producing a limited effect on the shape of the bandhead. 
In comparison, the consideration of non-LTE effects and spectral fitting of multiple bandhead 
regions provide powerful additional diagnostics of the detailed physical properties of the inner 
disk. Following a brief discussion of the evolutionary status of WL 16 (§3.4), we investigate 
these additional diagnostics. In § 4, we construct a non-LTE model and in §5 we apply the 
model to v=2-0 and v=5-3 observations of another YSO, 1548C27, in order to demonstrate the 
constraints that these additional considerations allow us to place on the physical properties 
of inner YSO disks. 

3.4- WL 16: On the Stellar Birthline or Main-sequence Star? 

In fitting the v=2-0 bandhead emission, we find that we can model the stellar component of 
WL 16 as either a G-star on the birthline or an A-star near the main-sequence. The presence of 
aromatic hydrocarbon features in the spectrum of WL 16 and the fit to the 2.3//m continuum 
provide other potential discriminants between these two possibilities. Of the various aromatic 
hydrocarbon features observed in the spectrum of WL 16, the 3.3/im feature is the strongest 
indicator of an early spectral type because it is excited by photons shortward of 3000A (Schutte 
et al. 1993). We can compare the efficiency required to excite the observed 3.3/im feature 
with UV radiation from either stellar spectral type with the efficiencies measured for other 
astronomical sources. For the models considered above, the fraction of the stellar (blackbody) 
luminosity that is emitted between 912A and 3000A, is 1% for T* = 5600K (Model L14) and 
6% for T* = 8700K (Model L15). For d = 160 pc, the strength of the 3.3/xm feature in WL 16 
is ~ 2 X 10 -3 Lq (cf. Tokunaga et al. 1991) where we have assumed an extinction correction 
at 3.3/im appropriate for Ak = 2. This luminosity is 0.9% of the available UV radiation for 
a 22Lq star with T* = 5600K and 0.14% for T* = 8700K. Although the latter value is more 
typical of those measured for Herbig AeBe stars, the former is not unreasonable given the 
large range in the measured values (Brooke et al. 1993). Thus, the strength of the 3.3/im 
feature does not strongly imply a spectral type earlier than G. 

While the continuum level is fit fairly well with a G-star on the birthline and an optically thick 
inner disk, if WL 16 is instead an A-star near the main-sequence, we require an additional 
source of 2.3//m continuum. The very small grains (VSGs) that are the larger counterparts of 
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aromatic hydrocarbons are likely to accompany these macro-molecules in the near environment 
of WL 16 and may contribute to this continuum. Natta et al. (1993) have shown that VSGs 
surrounding Herbig AeBe stars can reprocess stellar radiation into the near- and mid-IR 
(~ 2 — 20//m), significantly enhancing the strength of the continuum at 2.3/xm. Compared to 
the aromatic hydrocarbons responsible for the 3.3/xm feature, VSGs are excited by a much 
broader spectrum of photons. Thus they may significantly enhance the 2.3/xm continuum of 
WL 16 over a range of stellar spectral types and may make up the deficit of 2.3/xm continuum 
noted above for Model L15. 

In an effort to better determine the spectral type of WL 16, we obtained a high-resolution 
CSHELL spectrum centered on the 1.711/im Mg line. For stars of spectral type G and later, 
this line is among the strongest stellar absorption features in the H band. While WL 16 is 
substantially fainter at this wavelength (H = 10.5) than in the K band (A' = 7.8) the de- 
tectability of stellar lines should be better since the infrared excess in pre-main-sequence stars 
generally decreases with decreasing wavelength. The total integration time was 70 minutes 
through a 2" slit, giving a spectral resolution of ~ 23kms _1 and a signal-to-noise ratio of 70. 
No absorption features were present. The limiting factor on a line detection is the residual 
fringing in the spectrum which has a peak-to-peak amplitude of 4% and a period of about 
360 km s _1 . 

The non-detection of the Mg line could be due to any of the following reasons: the infrared 
continuum excess significantly dilutes the stellar spectrum; the stellar rotational broadening 
of the line is large; the star is hot enough that the line is weak or absent; or some combination 
of these. In order to investigate the constraints that can be placed on the spectral type of 
WL 16 from the if-band data, we calculated synthetic stellar spectra for the 1.711 /im Mg 
line. We used an updated version of the spectrum synthesis program moog (Sneden 1973) and 
the Kurucz grid of model atmospheres (Kurucz 1993). The (//-value and damping constants 
for the Mg line were determined by fitting an observed solar intensity spectrum (Livingston 
& Wallace 1991). 

In the G star model for WL 16 described in §3.2 (Model L14) ~ 70% of the 2.3/xm continuum 
comes from the star; the star should contribute an even larger fraction of the continuum at 
1.7/xm. For the assumed stellar parameters (T e ff = 5500K and logg = 3.5), the Mg line is 
predicted to be 35% deep at our spectral resolution in the absence of rotational broadening 
and other sources of continuum emission. In order for rotational broadening alone to account 
for the non-detection of the Mg line, we require vsini > 200kms _1 which is a significant 
fraction of the breakup velocity for such a star (310kms _1 ). A larger 1.7 /im continuum 
excess reduces the required rotational velocity. For example, with a more moderate vsini of 
50kms _1 we we require an excess flux of at least 3 times the stellar flux. We conclude that 
a late-type star cannot provide the major contribution to the near-infrared continuum unless 
the star is an extremely rapid rotator. 

The Mg line strength decreases with increasing effective temperature. For the temperature of 
our A star model (Model L15, T e ff = 8500K) the observed depth is predicted to be 17%. In 
this model, the A star contributes only 23% of the 2.3/xm continuum flux. Assuming that the 
star contributes ~ 35% at 1.7/im , a rather low vsini of 25kms _1 would be sufficient to hide 
the line in our data. Thus an A star, with the majority of the continuum coming from excess 
emission, would be consistent with the non-detection of the Mg line. 
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4. Non-LTE Model and Multi-wavelength Fits 

4-1. Level Population Calculation 

We calculate the population of the rovibrational levels using a ten-level CO molecule to rep- 
resent the vibrational levels and assuming that the rotational levels within each vibrational 
level are thermally populated. This approximation is based on the relatively large rate co- 
efficients for the collisional de-excitation of pure rotational transitions (Green & Thaddeus 
1976, Flower & Launay 1985, Schinke et al. 1985). In the temperature range of interest, the 
collisional de-excitation rate for both H and H2, k(v, J — >■ v, J — A J) ~ 5 X 10 -10 cm 3 s _1 for 
AJ = 1,2, . . ., is much larger than the corresponding rate coefficients for collisional changes 
in the vibrational (as well as rotational) quantum numbers. This approximation is valid even 
in the important case in which collisions with atomic hydrogen dominate the collisional rates 
of the rovibrational transitions, where k(v, J — > v — 1, J — A J) ~ 10 -11 cm 3 s _1 (cf. Appendix 
B). 3 The assumption of rotational equilibrium greatly simplifies the determination of the level 
populations and allows us to include in a simple way the effects of radiative trapping on the 
level populations. In this case, the vibrational level populations can be determined by the 
solution of a coupled set of effective vibrational level population equations which we derive in 
the following way. 

The large oscillator strengths and collision cross-sections of the Av = 1 transitions suggest 
that the population of a given vibrational level is primarily determined by radiative and 
collisional interactions with adjacent vibrational levels. Therefore, we begin by assuming that 
each vibrational level v is directly coupled only to the vibrational level above it, v' = v + 1. 
As discussed by Scoville et al. (1980), in this approximation successively higher vibrationally 
excited states are populated by climbing the vibrational ladder "one rung at a time". In 
steady state, the number of transitions into level v from level v' is balanced by the number of 
transitions into level v' from level v : 

Ju + C u i) = Y / n £(Bi u Ju + Ciu), (4) 

u,£ u,£ 

where we have adopted the notation £ = (u, J) and u = {v\ J') and the sum is over all transi- 
tions permitted between the upper and lower vibrational levels (e.g., the R and P branches). 
In this equation, n u = n(v\ J') and m = n(v } J) are the number densities of the upper and 
lower levels; A u i } B u i } and Bi u are the Einstein A's and B's of the transitions connecting 
the upper and lower levels; Ci u and C u i are the upward and downward collisional rates for 
these transitions assuming a sum over all collision partners (H, H2, e, etc.). We have also 
denoted the mean intensity of the radiation field weighted by the line profile function <f> v as 
J v = J J v cj) v dv. 



3 There is, in principle, an important difference between the rotational equilibrium of ground 
and excited vibrational levels stemming from the Einstein A- values for the two cases. Since 
the A-value is proportional to the cube of the transition frequency, the critical density for an 
excited vibrational level can be very large, typically n CI i t (v > 0, J) ~ 10 10 cm -3 . Nevertheless, 
the large difference between vibrational and rotational collisional rates means that the total 
vibrational populations are more readily disequilibrated, and we choose to focus on this aspect 
of a non-LTE calculation. 
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Consistent with our assumption of vertical homogeneity (e.g., of the level populations) in 

the line layer, we calculate J v at an "average" height in the slab, located where the vertical 

optical depth is half of the total. In addition, we ignore the effect of overlapping lines in the 

determination of the mean intensity. While the issue of overlapping lines is important for the 

calculation of the intensity emergent from the slab, we show below that the line overlap does 

not significantly affect the level populations. Thus, in the calculation of the level populations 

we treat each line as isolated and assume that the mean intensity which excites the line is a 

u i 

function of the optical depth for that line alone (i.e., J v = J v ). 
With these assumptions, J v for a given line is (cf. eq. I) 4 

J„ = ±- J ^ v dv J I u dtt = B„(T C ) ^ J ^ x dx e-^+'^dfi 

where all physical quantities such as B V (T C ) } S V} and ~f v are assumed to be independent of // 
because the disk is geometrically thin. The integration limits reflect the difference between 
the background radiation field which originates from the 2tt steradians subtended by the 
continuum layer and the diffuse field which originates from all Att steradians. The optical 
depth 

T v = J T v dv (5) 

is an average over the line profile where Av is the characteristic width of the line profile 
centered on the line frequency vq. 

We introduce a dimensionless line profile function cf> x = cf> v Av such that t v = ~T v (f> x and x = 
{y — v§)j Av . Following Mihalas' treatment of the escape probability in a static plane parallel 
atmosphere with vertically homogeneous properties (Mihalas 1978), we interpret quantities of 
the form 



/OO fil ft OO 

<P x dx / e- T +'l*dii= \ E 2 (T<p x )<p x dx (6) 
- oo J J — oo 

as escape probabilities. The line profile function in the outer integral expresses the probability 
that at a vertical optical depth r a photon is emitted x = (v — v§)j ' Av line widths from line 
center. Such a photon has a probability E2(r(f> x ) of escaping the layer in any direction fx. Our 
expression for /3(r) differs from Mihalas' expression for Peij) by a factor of 2 because we have 
assumed that the photon can escape out of both the top and bottom of the atmosphere and 
that the optical depth in either direction is r = T v /2 on average. While the photon cannot 
really "escape" by propagating into the continuum layer below, we imagine that photons 
reaching the continuum layer are destroyed and reemitted as continuum radiation. With this 
identification, the mean intensity weighted by the line profile is 

l v = B v (T c )P(t v /2)/2 + S v [l- P(t v /2)] . (7) 



4 



In the interest of notational economy, we have dropped the superscripts ui on quantities such 
as J v , lv , *5*j/, 7~v and j3. 
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Note that the limits of integration in dimensionless frequency x in equation (6) assume an 
isolated line. For the CO fundamental transitions, the P-transitions are separated by 4- 
8 cm -1 and the R-transitions are more closely spaced (0-6 cm -1 ) since the transitions go 
through a bandhead. If we conservatively estimate the average spacing between neighboring 
lines to be Av = 1cm -1 , truncating the integral at x t = Atj/Stj (~ 50 for a linewidth 8v of 
2.5 km s -1 ) reduces the escape probability by less than 50% relative to the untruncated case. 
So overlapping lines do not have a dominant effect on the relevant optical depths and escape 
probabilities. 

Using the definition of the line source function to rewrite equation (7) for a given line as 



— ul 
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e hv/kT c 



where f3 u i = /3(r" £ /2) and ni, n u , gi, and g u are the number densities and statistical weights 
of CO in the lower and upper levels, equation (4) becomes 
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If the rotational levels are thermally populated at the gas kinetic temperature Tt 
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yj_ _ (2 J + 1) J( J+ i)B/ T/ = f 



(9) 



where n v j = n(v } J), n v = ^2j n vJ} B is the rotational constant, and Z ~ Tt/B is the 
partition function for the rotational levels, the radiative terms in equation (8) can be further 
approximated as 

(n v /n v i - 1) /2 



e hu/kT c 
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With this approximation, the assumption of thermally populated rotational levels, and a little 
algebra, equation (8) becomes 
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(10) 



where 



Pv'v = 2 X/ f J ' @ ul ' 
ul 



ul 

C'vv 1 = X^ f J Cul ' 



(11) 



The terms A*, v , f3*, v7 C*, v7 and C* v , } can be regarded as effective Einstein A's, escape prob- 
abilities, and collisional terms that relate the vibrational levels. Note that we have made the 
additional approximation of separately averaging the A-values and escape probabilities since 
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the A-values usually do not vary substantially over the range of J in which the escape proba- 
bilities are substantial. The factor of 1/2 in the definition of f3*, v arises from the consideration 
of the R and P transitions out of the upper levels. 

By detailed balance the upward and downward collision rates satisfy C j> j = (nj/nji )lte Cjji. 
The assumption of thermal rotational level populations now implies that 

r<* — (v'-v)e/T r <* 

where 6 is the vibrational constant. When collisions dominate, equations (10) and (12) yield 
the LTE result: n v > /n v = exp[— (v' — v)8/Ti]. When radiative transitions dominate, equation 
(10) implies n v i /n v = l/[2 exp(hiy/kT c ) — 1], which differs from the Boltzmann relation for 
excitation temperature T c because the underlying radiation field is, as noted earlier, only half 
a blackbody. 

The system of coupled, statistical-equilibrium equations for the vibrational level populations 
n v is linked to the radiative transfer part of the problem through the determination of the 
escape probabilities f3 u i = /3(r"^/2) which requires the self-consistent determination of the 
vibrational level populations that enter in the optical depths r" £ . This suggests that the system 
of equations can be solved iteratively given the initial assumption of LTE level populations. 
The rates for vibrational excitation due to collisions with H and H2 are based on the rate 
coefficients discussed in Appendix B. The total collision rates, C*, v and C* v ,, are obtained 
by assuming chemical equilibrium abundances of H and H2. Although Einstein A-values for 
individual lines within a vibrational band can vary by factors of several, the average A-value 
between two vibrational levels A^, v varies weakly with temperature (see also Appendix A). 
Therefore, we adopt a temperature-independent rate of A^, v = A* /t) (T = 2000K). The set 
of equations (10) can easily be generalized to couple more than two neighboring vibrational 
levels together. It is usually not necessary to couple more than three levels together to solve 
the set of equations (10) to a fractional accuracy of < 10 -7 . 

4-2. Spectral Synthesis of the v=5-3 and v=2-0 Bandhead Regions 

The u=5-3 bandhead region of the spectrum provides valuable complementary information 
to that contained in the u=2-0 bandhead region. For example, the onset of non-LTE level 
populations is difficult to detect with measurements of the u=2-0 bandhead alone. This is 
evident from a comparison of the fits to the WL16 u=2-0 bandhead region shown in Figures 3b- 
c and 3e-f. The latter set of models fits, obtained with the non-LTE model described in the 
previous section, show that similar fits to the data can be obtained with the nearly the same 
parameters as in the LTE case (cf. Table 1). This is because the u=2-0 bandhead region 
consists of lines that are primarily sensitive to gas at temperatures larger than those at which 
departures from LTE become important. This is shown graphically in Figure 4 in terms of the 
relative optical depth of a line at the u=2-0 bandhead in LTE as a function of temperature. 
Because lines in the u=2-0 bandhead arise from relatively high-lying rotational levels, they 
are not well excited below 1500 K even in LTE. Thus, while the lines in the u=2-0 bandhead 
provide useful information about the kinematics of the gas in the inner disk, they probe only 
a limited range in temperature. 

In comparison, isolated lines from several other CO overtone transitions also fall in the u=5-3 
bandhead region (Fig. 5); these lines and the 5-3 bandhead together probe a much wider range 
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of excitation temperatures (500K<T<5000K) than those contributing to the v=2-0 bandhead. 
In Figure 4, the relative LTE optical depths fi u (2Ji + l)exp(—Et/kT) of representative lines 
from different overtone transitions in the v=5-3 bandhead region are compared with that of 
a line at the v=2-0 bandhead. The strongest lines in the v=5-3 bandhead region are the 
2-0 and 4-2 lines and the 5-3 bandhead itself. The 2-0 P transitions are the primary low 
temperature constraint. Because the 4-2 and 2-0 lines do not strongly overlap, it is possible 
to establish their relative strengths. The relatively large separation in wavelength between 
these lines also makes it possible to measure the line profiles of individual lines for sources 
with projected rotational velocities 100 kms -1 . In the next section, we present our high 
spectral resolution observations and modeling of the v=2-0 and v=5-3 bandhead emission 
from such a source, the probable Herbig AeBe star 1548C27. 

5. / =5-3 and / =2-0 Bandhead Spectroscopy of 1548C27 

5.1. What is 1548C27? 

Discovered as a cometary nebula (Craine et al. 1981), the YSO 1548C27 has an associated 
jet which is aligned with the nebular axis (Mundt et al. 1984). Energetic outflow activity is 
also indicated by the optical spectrum of the central object which is characterized by P Cygni 
Ha, H/3 and broad blueshifted Nal D absorption (Mundt et al. 1987). In the near-infrared, 
1548C27 has very prominent CO overtone and Br 7 emission (Carr 1989, 1990). The [/-band 
- 100 /im spectral energy distribution indicates a system luminosity of L ~ 128(d/kpc) 2 Lq 
(Vilchez et al. 1989) with the spectral shape indicating Herbig AeBe Group II membership 
(cf. Hillenbrand et al. 1992). 

The distance to 1548C27 is controversial. It is located in front of the HII region S86 and in 
the same region of the sky as NGC 6823, both of which are located at a distance of ~ 2.5 - 
3 kpc. Arguments based on the spectral energy distribution shortward of 3 /im, suggest that 
1548C27 is at a distance of d = 0.91 kpc while CO ( J=2-l) and CO ( J=3-2) data have been 
used to derive a kinematic distance of 2.4 kpc (Dent & Aspin 1992). 

5.2. CO Overtone Spectroscopy 

Spectra of 1548C27 in the v=2-0 and v=5-3 bandhead regions were obtained on 1993 May 
15 at the NASA Infrared Telescope Facility using the facility cryogenic echelle spectrograph 
(CSHELL; Tokunaga et al. 1990, Greene et al. 1993) with a NICMOS 3 256x256 HgCdTe 
array and a 1" slit width that provided a resolution of 14 kms -1 . The data were acquired 
and reduced as described in Carr et al. (1993). The resulting spectra for the v=2-0 and 
v=5-3 bandhead regions are shown in Figure 5. As in the case of WL 16, the 2-0 bandhead 
of 1548C27 also shows the characteristic shape of bandhead emission from a rotating disk. 
Although the velocity extent of the bandhead emission (100 — 120 kms -1 ) is less than in the 
case of WL 16, the characteristic blue wing, shoulder and red peak are clearly present. 

In the v=5-3 bandhead region, the feature that peaks at 2.3837/im is a complex blend of 
the 5-3 bandhead, located at 2.3829/im, and neighboring 4-2 R and 2-0 P lines. Although 
the kinematics of the emitting gas are difficult to deduce from the shape of this feature, the 
double-peaked shape of the features at 2.3795 and 2.3810/im strongly support the rotating 
disk interpretation for the emission. The feature at 2.3795/xm is essentially a 4-2 R line and 
the feature at 2.3810/xm is a blend of a 4-2 R line and a 2-0 P line. Based on the velocity 
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extent of the v=2-0 bandhead emission, these features overlap only slightly in their respective 
red and blue wings and are a reasonably good representation of an isolated line profile. 

5.3. Modeling of the v=2-0 and v=5-3 spectra 

To demonstrate the constraints that these observations place on the inner disk properties of 
the system, we assume that the distance to 1548C27 is 2.4kpc (cf. Dent & Aspin 1992) and, 
therefore, the system luminosity is 740^0. Studies of the evolution of intermediate mass pre- 
main sequence stars restrict the stellar mass of such a system to ^ 5Mq (Palla & Stahler 
1993). For this mass, the velocity extent of the 2-0 bandhead implies that the system must 
be viewed at relatively low inclination. For a A'-band extinction of Ak = 0.62 mag (Vilchez 
et al. 1989), and i = 30°, we can roughly match the velocity extent of the bandhead emission 
with M* = 4Mq and a CO dissociation radius at ~ 2Oi?0. Assuming that the continuum 
layer has the usual radial temperature dependence T c = T c o(r/ro) _3/ ' 4 , we can fit the 2.3 //m 
continuum level and the system luminosity with i?* = ro = 5.9i?0 and T c o = 9800K. These 
parameters correspond to a stellar luminosity of = 200 Lq, an age of 1.5 X 10 5 yr from the 
birthline (Palla & Stahler 1993), and a disk luminosity of Lr> = 47ri?* 2 crT ( 5 = 280 Lq. We 
assume that a boundary layer/hot spot contributes another 260 Lq so that the total system 
luminosity is 740 Lq. 

The adopted stellar parameters place 1548C27 among the Herbig AeBe stars. In order to place 
a lower limit on the disk accretion rate, we conservatively assume that 25% of the non-disk 
luminosity (460 Lq) is reprocessed by the disk; the disk then has an accretion luminosity of 
Lr> acc = (280 — 115) = 165 L©. 5 This value of Lr> corresponds to a disk accretion rate 
of Mr> ~ 1 X lO -5 M0yr -1 (eq. 3) which is consistent with the evidence for energetic mass 
outflow from 1548C27 and with the large disk accretion rates characteristic of Herbig AeBe 
stars. 6 

Applying the non-LTE model with these parameters, we obtain a good fit to the detailed 
properties of both spectral regions with Lorentzian line profiles and the parameters shown 
in Table 2. In each case, the synthesized spectrum for the 2-0 bandhead region is scaled 
by ~ 95% to match the continuum in both spectral regions (Fig. 6). In the 5-3 bandhead 
region, the true continuum is 0.246 Jy whereas the apparent continuum of 0.26 Jy (e.g., as 
measured at 2.382/im ) is formed by the overlap of rotationally broadened lines. The "wiggles" 
in the emission redward of the peak in the 2-0 bandhead region arise from the "beating" of the 
(double-peaked) rotational broadening function against the rest wavelength distribution of the 
overtone lines. For systems with reduced projected rotational motion (e.g., 1548C27 and SVS- 
13), this effect appears close to the 2-0 bandhead where the spacing of lines is comparable 

5 This probably underestimates the actual accretion luminosity because the disk albedo is prob- 
ably nonzero and because the reprocessed radiation should not be counted separately from 
that of the primary source in computing the total system luminosity. Note, however, that the 
deduction of L sys from observations depends on the viewing angle. 

6 Because the CO dissociation radius is larger than the stellar radius and the stellar contribution 
to the 2.3/im continuum is negligible, it is possible to find other stellar parameters which do 
not substantially alter the fit but correspond to different system luminosities. Note also that 
if the distance to 1548C27 is 0.9 kpc (Vilchez et al. 1989), the system luminosity is reduced 
by a factor of 7 and the disk accretion rate could be as low as 1O -6 M0 yr -1 . 
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to the separation of the horns of the rotational broadening function. The high resolution 
overtone spectroscopy by Carr & Tokunaga (1992) and Chandler et al. (1993, 1995) show the 
same effect. 

As shown in Figure 6, we are able to fit the 2-0 shoulder, peak, blue wing, and slope of the 
emission redward of the peak reasonably well with 7 in the range 1.0 kms -1 (Model Gl_10) 
to 5.0 kms -1 (Model G5_9). The model "wiggles" would have reduced amplitudes if we 
convolved the synthetic spectra with the instrumental profile, but they would still appear too 
prominent when compared with the observations. This indicates perhaps a smoother inner 
cutoff to the CO abundance than given by the chemical and mechanical equilibria calculations 
of this paper. Apart from this problem, the relative heights of the emission at 2.3800 /im and 
2.3815 /xm are well fit by both models, as is the shape of the 5-3 bandhead at 2.3829 /im. 

With 7 = lkms -1 , there is a relatively narrow range of Tt and at within which acceptable 
fits can be obtained. If the outer radius to the emission is roughly held fixed, much flatter at 
and steeper Tt distributions would steepen the slope redward of the 2-0 bandhead and weaken 
the 2-0 wing relative to the fit in Model Gl_10. This is because both a flatter at and a steeper 
Tt reduce the emission from small radii (the latter by dissociation), weakening the blue wing 
and steepening the slope redward of the bandhead. Although equally good fits to the 2-0 
bandhead are possible for somewhat flatter Tt and at distributions, these parameters would 
also overemphasize the strength of the 2-0 P lines in the u=5-3 bandhead region. Temperature 
distributions as flat as Tt ~ r -0 ' 3 can be ruled out by the shape of the 2-0 bandhead alone. In 
this case, the flatter Tt distribution reduces the emission from small radii by reducing the line 
source function, steepening the slope redward of the bandhead and weakening the blue wing. 
Intrinsically broader lines that are less optically thick allow acceptable fits within a larger range 
of Tt and a t distributions. Model G5_9 with T t = 9200(r / i?*) -0 - 55 and a t = 620(r/i?*)- 2 
also produces a good fit. Column density distributions as flat as r _1 can fit both bandhead 
regions reasonably well with Tt ~ r -0 ' 6 . 

5.4- Discussion of the Non-LTE Effect 

In these models, the strength of the feature at 2.381 /im and the shape of the bandhead at 
2.384 /im are sensitive to the departure of the level populations from LTE beyond the radius 
at which Tt drops below 2000K. This effective outer radius to the emission occurs somewhat 
beyond the radius at which H associates to form H2 , since H2 has a cross-section for collisionally 
exciting CO that is smaller than that of H by almost two orders of magnitude (see Appendix 
B). When non-LTE effects set in, the upper vibrational levels are depopulated first, followed 
by lower vibrational levels. This sequence produces effective outer radii that are larger for the 
lower overtone transitions. Radiative trapping plays an equally important role in determining 
the radii at which these effects occur because the line layer column densities required to fit 
the overtone emission correspond to very large optical depths in the fundamental transitions. 
The low escape probabilities associated with these optical depths result in significant radiative 
trapping of fundamental transition photons (cf. Scoville et al. 1980). 

At these large optical depths, photons can only escape in the wings of the line, so the degree 
of radiative trapping depends significantly on the assumed behavior of the wings of the line 
profile function. For example, the escape probability for gaussian lines depends asymptotically 
on the line optical depth as ~ r" 1 (logr^) whereas the escape probability for lines with 
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Lorentzian or Voigt profiles depends asymptotically on the line optical depth as ~ t v (see 
Mihalas 1978). The shape of the line profile function also affects which transitions contribute 
most significantly to the escape of fundamental photons. 

At the column densities required to fit the overtone emission, escapes from line profiles with 
weak line wings, such as gaussian line profiles, are difficult because they must be made far out 
in the wings. As a result, only the high- J i transitions have any substantial escape probability 
due to the reduced population of these levels and their consequently lower optical depths. At 
the same time, since the population of these levels is small, their contribution to the effective 
escape probability (eq. 10) is also small because the weighting factor (fj>) is also reduced by 
the level population. This radiative trapping results in the extension of LTE level populations 
to radii significantly beyond the radius at which H associates to form H2 . 

In comparison, the probability of escape from lines with Lorentzian profiles is large enough 
that the transitions from relatively well populated J 1 levels contribute significantly, resulting 
in substantially larger effective escape probabilities. In the fit shown in Figure 6a-b, which 
has a Lorentzian line profile with a width of 7 = 1.0 kms -1 , the vibrational level populations 
reach half their LTE value at ~ 1800K (v=5) to ~ 1400K (v=l). Correspondingly, the strong 
transitions in the v=5-3 bandhead region have source functions that depart from LTE at 
~ 1800K (u=5-3), ~ 1750K (u=4-2), and ~ 1550K (u=2-0), introducing effective outer radii 
to the CO emission at these temperatures. In this particular case, this departure from LTE 
only affects the spectral fit through its impact on the strength of the 2-0 P lines in the v=5-3 
bandhead region. This is because only the 2-0 P lines still have a significant optical depth at 
the temperature where the departure from LTE is significant: at ~ 1800K the optical depth 
of a representative 5-3 line is r(5-3) ~ 0.2, at ~ 1750K r(4-2) ~ 3, and at ~ 1550K r(2-0) ~ 
11. 

To illustrate the magnitude of the non-LTE effect on these lines, in Figure 7a-b we also 
show the spectral fit obtained by the (artificial) extension of LTE level populations to 600K 
(r = 100i?*). The synthesized v=2-0 bandhead is somewhat stronger at the peak and the 
"wiggles" are filled in by emission from larger radii. The overall shape of the v=2-0 bandhead 
is altered only slightly, because the lines in the v=2-0 bandhead have optical depth r ^ 1 at 
the radii at which the non-LTE effect would have set in. In contrast, the overall shape of the 
spectrum in the v=5-3 bandhead region is altered significantly by the prominence of the 2-0 
P lines which have r > 1 to r > 40i?* are consequently enhanced due to emission from larger 
radii. This sensitivity illustrates the utility of the 2-0 P lines as diagnostics of low temperature 
gas. To illustrate the effectiveness of this non-LTE effect in truncating the emission beyond a 
given radius, we show in Figure 7c-d the spectral fit obtained with non-LTE level populations, 
the parameters of Model Gl_10, and an abrupt truncation of the emission at the radius at 
which T = 1500K. 

6. Discussion and Summary 

We have successfully modeled the CO overtone bandhead regions of the observed spectra of 
the YSOs WL 16 and 1548C27 in terms of emission from Keplerian disks. Our observations 
provide some of the most compelling evidence to date for the existence of disks around young 
stars. Through the detailed study of this emission, we are able to place useful constraints 
on the detailed physical properties of the inner disk region. These results also indicate the 
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possibility of dynamical mass estimates of young stars through the study of the emission 
from their surrounding disks. In the case of WL 16, we are able to place a useful constraint 
on the stellar mass due to the relatively high inclination of the system and the availability 
of additional observational constraints on the distance and system luminosity. In the case 
of 1548C27, observations of the v=2-0 and v=5-3 bandhead regions provide a particularly 
potent combination of diagnostics which reveal the onset of non-LTE level populations. In 
our non-LTE model for the emission from this source, the CO emission arises from gas at 
temperatures between ~ 5000 — 1500K. The inner radius is due to the dissociation of CO; the 
outer radius is the result of the departure of the vibrational level populations from LTE due 
to radiative escapes from suprathermally broadened lines. The need for this outer radius is 
clearly indicated by the strength of the 2-0 P lines in the v=5-3 bandhead region. 

Thermal Doppler line profiles appear unable to reproduce either the shape of the v=2-0 and 
v=5-3 bandhead regions of 1548C27 (§5.3) or the strength of the v=2-0 bandhead emission 
from WL 16 given our adopted extinction to the source (§3.2). In contrast, Lorentzian line 
profiles with widths of 7 = 0.5 — 5 kms -1 produce excellent fits to the overtone bandhead 
regions of both sources. These considerations argue against the presence of dust in the tem- 
perature inversion region of the disk atmosphere as the sole explanation for the effective outer 
radius to the overtone emission. In the absence of turbulence, dust is believed to settle very 
quickly from the surface layers of a circumstellar disk (Weidenschilling & Cuzzi 1993). In 
the presence of turbulence, which may mix dust into the surface layers, the non-LTE effect 
discussed in the text will come into play independent of any radiative transfer role for the 
dust. 

The relatively large inner disk temperature inversion and the steep temperature and column 
density distributions required to fit the overtone emission restrict the possible physical expla- 
nations for the origin of the inversion. These include irradiation of the disk surface by the 
star and funnel flow accretion shock; turbulent heating in a wind-disk interface; and chromo- 
spheric heating that results from the vertical propagation of MDH waves into the upper disk 
atmosphere. Calvet et al. (1991) have used an approximate irradiation model to explore the 
ability of stars to heat their circumstellar disk surfaces and whether the CO overtone band- 
heads consequently appear in emission or absorption as a function of stellar spectral type and 
disk accretion rate. However, in terms of the detailed properties of the inner disk region, the 
model predicts flatter temperature and column density distributions as well as smaller tem- 
perature inversions than are implied by our observations of WL 16 and 1548C27. The steeper 
temperature and column density distributions required to fit the observations could result, 
instead, from the interaction of a stellar wind with the circumstellar disk: as the wind blows 
over the disk, the disk atmosphere is stirred up and heated by turbulent dissipation, forming 
an aeolosphere (cf. Carr et al. 1993). Suprathermal widths may also characterize lines formed 
in disk chromospheres. If accretion in the inner disk is the result of the Balbus-Hawley insta- 
bility (Balbus & Hawley 1991), the magnetic fields that participate in this instability may also 
propagate MHD waves into the upper disk atmosphere, producing a turbulent temperature 
inversion layer. 

If the turbulent line broadening resulting from either of the latter two possibilities is similar 
in nature to that observed in spectral line profiles of turbulent interstellar gas (Falgarone & 
Phillips, 1990), the Lorentzian line profiles adopted in our synthetic spectra may reasonably 
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represent the expected non-gaussian line wings. We have examined and rejected the alterna- 
tive possibility that the enhanced line widths are due to pressure broadening. According to 
standard estimates for this broadening mechanism, the gas pressure in the model layers are 
too low by several orders of magnitude to be responsible for the line widths. We conclude 
that the line-emitting region in these sources is likely characterized by velocity widths much 
larger than those appropriate for a quiescent layer heated by the external radiation field of 
the star. Since external irradiation will tend to stabilize disk atmospheres against convection, 
we suggest that these turbulent motions arise from perturbing influences such as stellar winds 
or vertical MHD wave propagation. 

The stellar properties and disk mass accretion rates implied by our observations also bear 
on the origin of the inner disk temperature inversion. In the case of WL16, the implied 
disk accretion rate is low, Mo ^2 X lO -7 M0yr -1 , independent of whether we model the 
stellar component as a G-star close to the stellar birthline, or as an A-star close to the main 
sequence. Given the low disk accretion rate, heating due to irradiation by the relatively 
early type star may then provide a simple explanation for the temperature inversion in the 
upper atmosphere of the surrounding disk. The low disk accretion rate is also consistent with 
the lack of evidence for energetic outflow activity from this source; at such low wind mass 
loss rates, aeolospheric or chromospheric effects may be confined to a modest amount of line 
broadening which contributes to the strength of the v=2-0 bandhead emission. However, the 
steep temperature and column density distributions derived for the line-emitting region may 
indicate that these effects are also responsible for additional non-negligible heating. 

In contrast, our models for 1548C27 indicate a large disk accretion rate, Mo ~ 1 X 10 -5 Mq yr _1 .| 
Although this large mass accretion rate is consistent with both evidence for energetic mass 
outflow from 1548C27 and the disk accretion rates of other young stars with stellar parameters 
similar to those deduced for this object (Herbig AeBe stars), it also places 1548C27 in a pa- 
rameter regime in which the models of Calvet et al. (1991) would nominally predict that the 
CO lines should be in absorption rather than in emission. This discrepancy may be reduced 
by considering the impact of the funnel flows that are likely to accompany these large disk 
accretion rates if intermediate mass stars possess magnetospheres (cf. Najita et al. 1995): the 
UV emission from the accretion hot spots at the base of the funnel flows is an energetic source 
of radiative heating. The possibility of a large contribution to the temperature inversion of the 
inner disk from funnel flow hot spots introduces the possibility of another source of energetic 
heating for the inner disk. As described by Shu et al. (1994), stellar magnetic fields strong 
enough to truncate the disk and to give rise to a funnel flow are also strong enough also to drive 
a magnetocentrifugal X-wind from the inner disk region. An aeolosphere may be created as 
the lowermost streamlines of this induced X-wind sweep over the inner disk region, producing 
the suprathermally broadened lines and steep temperature and column density distributions 
implied by our modeling of the overtone emission from this source. Similar line broadening 
and inner disk heating may arise from a disk chromosphere. 

An improved understanding of the origin of the inner disk temperature inversion requires a 
detailed theoretical study of the temperatures and column densities that characterize aeolo- 
spheres and disk chromospheres. In addition to indicating the degree to which aeolospheres 
and chromospheres can heat the upper disk atmosphere, these studies may also indicate the 
observational signatures of aeolospheres and chromospheres. When combined with observa- 
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tions of the type presented in this paper, these distinguishing characteristics could be used 
to address other astrophysical issues. For example, the observational identification of aeolo- 
spheres as the origin of the disk temperature inversion would favor a stellar origin for the wind 
over a disk origin, whereas the identification of chromospheres as the origin of the temperature 
inversion would support the Balbus-Hawley instability as the physical process that gives rise 
to accretion in disks. 
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AST-9114940, and NASA grant NAGW-2278 (J. S. C. and A. T. T.). This work was carried 
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Appendix A. Radiative Rates 



Approximate A-values 

The Einstein A-values for individual fundamental transitions from an upper level (u',J') to a 
lower level (v y J) can be obtained from 

6471"^ I Tfl I 

A v ,j,, v j = -^^(v'J'H^lvJfjjr^j (Al) 

given the transition frequency v and the electric dipole moment function fi(x) as a function 
of the reduced internuclear distance x. The quantity m = [J'(J' + 1) — J(J + l)]/2. For the 
purpose of calculating the vibrational level populations, we approximate the expectation of 
fj.{x), M% = (v' J'\[i(x)\v J), with the rotationless transition dipole moment M% (0): 

A -^=^ M "'^M-i (A2) 

and use transition frequencies v calculated using energy levels from Mantz et al. (1975) and 
values of M% (0) from Goorvitch & Chackerian (1994). 

Values of A*, v (cf. eq. 11) can be calculated by adding the rates for the R- and P-transitions 
associated with each upper level, weighting by the population of the upper level, and summing 
over all upper levels: 

v' j',vJ'-i + A v iji, v ji +1 ). 

J' 

For the purpose of calculating A*, v , we have approximated the rotational level populations as 

A ' J = ^^ eXp( " J(J + 1) ^ /T) 

where the rotational constants B v are from Mantz et al. (1975). Since A*, v varies weakly with 
temperature over the range T = 200 — 8000K, we further approximate A*, v with its value at 
2000K. 

Approximate Escape Probabilities 

The escape probability of fundamental photons averaged over rotational levels (eq. 11) are 
approximated in the following way. The Einstein B-values for the fundamental transitions are 
related to the A-values by 

B - ° 2 A 2J, + 1 

BvJ > v,JI -^ Av,ji > vJ JjTT- 

With the approximate form for A v >j^ v j given by equation (A2), 

| Tfl | 

B v j jV 'j' = b v i v 2 j _|_ (A3a) 
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where 

b v , v = (327r 4 /3/i 2 c)[M;'(0)] 2 . (A36) 

The sum of the B-values for the R- and P-transitions from the same lower level is independent 
of J 

B v j :V i j+i + B v j :V ij-i = b v i v . 
Therefore, b v > v is the effective B-value relating the vibrational levels v and v' 

B* vv , = b v > v . (A3c) 

The corresponding optical depths are obtained from these upward B-values (eq. A3) using 

vre 2 _ hv \m\ 
-Jlu — — & vv i 



m e c^ 4vr 2 J + 1 

for the absorption oscillator strength fi u in the usual expression for the line optical depth. 
We also use the following simple approximation for the vibration-rotational energy levels in 
order to calculate transition frequencies and stimulated emission terms: 

E vJ = v8 + J(J + 1)B , 

where 6$ = 3084K and Bo = 2.768K. Accordingly, for the rotational level populations, we 
adopt 

fj = { if/B~) exp(_ J(J + 1)Bo/T) - (A4) 

Given these optical depths, the effective escape probability f3*, v (cf. eq. 11) is determined by 
averaging the escape probabilities for the R- and P-transitions associated with each upper 
level, weighting by the population of the upper level, and summing over all upper levels: 



Pv'v = ~ ^2 fj'(Pj',J'-i +/3j',J' + i) 



2 

J' 



where = (3(tJ ,j 1 ), (3ji,ji+i = /3(t/ ,j +1 ), and the rotational level populations are 

as approximated in equation (A4) 



Appendix B. Collisional Excitation of CO Vibrational Levels 

The excitation of the rotation- vibration bands of CO has been discussed earlier in a variety 
of astrophysical contexts, e.g., Thompson (1973; late-type stars); Scoville, Krotkov, & Wang 
(1980; young stellar objects); and Ayres & Wiedemann (1986; the sun). Our discussion focuses 
mainly on collisional excitation by atomic hydrogen, where there are real differences between 
authors, in contrast to collisions with electrons and hydrogen molecules. We emphasize new 
experimental data on the role of Av > 1 transitions and the results of Elitzur (1983) based 
on surprisal theory. 
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Ayres & Wiedemann (1986) gave a critical discussion of existing experiments and theory 
on H + CO vibrational excitation. They used the Landau- Teller parameters determined in 
the shock-tube experiment of Glass & Kironde (1983, henceforth GK) and standard theory 
(Herzfeld & Litovitz 1959) to obtain the Av = 1 — de-excitation rate coefficient k(l — 0, H). 
GK found that the product of the post-shock pressure P and the measured relaxation time r 
is essentially constant in the temperature range of the experiment T = 1000 — 3000 K, 

Pt = 1.8 ±0.4 x 10" 8 atm-s. (Bl) 

Substituting this into the standard theory for the relaxation time 

k(l - 0) = k B T/Pr (l - e-* /T ) 1 (B2) 

we get a useful, direct representation of the GK data, 

&(l-0, H) = 7.57 x 10" 15 cm 3 sec" 1 T(l - e^ ^)" 1 . (B3) 

It must be emphasized that GK measured only one relaxation time and that they used the 
traditional method of analysis, based on the Landau- Teller theory, where the collision rate 
coefficients obey the rules for dipole matrix elements of a simple harmonic oscillator (the 
same as for radiation), i.e., there are no Av > 1 transitions. The GK result for H + CO, i.e., 
equation B3, is 100-20 times larger the rate coefficient they measured for H2 in the temperature 
range from 1000-3000 K. 

At the time of the Ayres & Wiedemann (1986) paper, there was little experimental information 
on the validity of the harmonic oscillator selection rules, which forbid Av > 2 transitions. 
Hooker & Millikan (1963) did study the problem for rare-gas and molecule collisions with 
CO with shock tube experiments at 2000 K. From the characteristics of the fundamental and 
first overtone time signals, they concluded that the CO population is largely built up by a 
sequence of Av = 1 collisions, i.e., by "climbing the vibrational ladder". They estimated 
that the Av = 2 collisions cannot occur as frequently as 1/10 as often as Av = 1 collisions. 
In a high-pressure experiment on vibrational relaxation of electronically excited CO in the 
A 1 ]! level, Fink & Gomes (1974) found that Av = 2,3 collisions were almost as important as 
Av = 1 collisions. However, it is likely that physical processes other than those operative for 
the ground electronic level are at work in this experiment. 

A new generation of molecular beam experiments is beginning to provide information on the 
dependence of the H + CO collisional excitation rate coefficients with selected final vibrational 
and rotational quantum numbers (e.g., Wight & Leone 1983, McBane et al. 1991). Although 
the experiments have only been done so far at energies higher than are relevant for astro- 
physical applications involving CO (1-3 eV), they do indicate that Av > 2 transitions are 
significantly smaller than Av = 1 transitions and that they decrease in importance as the en- 
ergy is decreased. This is in conflict with Elitzur (1983), who correctly criticized the complete 
ignoring of Av > 2 transitions, but proposed rate coefficients which we believe are too large, 
especially for transitions involving large changes in v. Elitzur 's work is based on a statistical 
theory known as "surprisal theory" (e.g., Levine & Bernstein 1976). We suspect that the data 
underlying his application of this method is inappropriate to the H + CO system, i.e., they 
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refer to different systems and situations, e.g., the one piece of evidence cited for CO involves 
rare gas atoms de-exciting an excited electronic level of CO (Fink & Gomes 1974) 

In this paper, we adopt the following rate coefficients: 

Electron Collisions - The rate coefficient for T < 5000 K that we derive from laboratory 
experiments (e.g., Morrison 1988) is 

k(l - 0,e)~1.24 x 10" 11 cm 3 s _1 T 1/2 . (B4) 

On this basis, we can ignore electronic collisions compared with neutrals collisions for x e < 
10 -3 , a restriction satisfied by the disk atmospheres and winds considered here. 

Molecular Hydrogen and Helium - We follow previous studies and use the results of shock-tube 
experiments expressed in Landau- Teller form, 

k(l - 0, X) = k B T e -(*T- 1/S -B) ^ _ e ~0/T^ - 1 atm _ g (B5) 

The parameters given by Milliken & White (1963) are: X = He, A = 99 atm-s and B = 20.4 
atm-s; X=H2, A = 68 atm-s and B = 19.1 atm-s. 

Atomic Hydrogen - In accord with the above discussion, we use equation B3. 

The available information suggests that Av > 2 transitions are relatively unimportant al- 
though quantitative data is not yet available for the temperature range of interest. We ignore 
these collisions and, again consistent with available experimental information, assume that 
all Av = 1 transitions have the same rates. The rate coefficients summarized above indicate 
that for the disk layer analyzed in this paper, only atomic and molecular hydrogen collisions 
are important. Accordingly, in the level population calculation described in §4, the total 
vibrational collisional excitation rate is 

nC;, v = n H £(1-0, H) + n H2 £(1-0, H 2 ). (B6) 

This rate depends on the ratio njj/^H 2 which is assumed to be the chemical equilibrium value. 
Given the A-values discussed above, the corresponding critical densities for the vibrational 
levels are rin ~ 10 12 — 10 13 cm -3 over the temperature range T = 2000 — 5000K. The critical 
densities increase toward both lower temperatures (due to reduced rate coefficients) and higher 
densities (which favor H2 over H). 
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Table 1. WL16 Model Parameters 



Model* 


R* 




Pc 




Pi 






<f>v 




(Rq) 


(K) 




(K) 




( g cm -2 


) 




L12 


5.0 


3000 


0.75 


6000 


0.75 


5000 


4.5 


Thermal Doppler 


L14 


5.0 


3000 


0.75 


6000 


0.78 


370 


3.8 


Lorentzian (7 = 0.5) 


L15 


2.08 


3000 


0.75 


5810 


0.78 


310 


3.7 


Lorentzian (7 = 0.5) 


U2 


5.0 


3000 


0.75 


5500 


0.75 


200 


2.8 


Lorentzian (7 = 0.5) 


V2 


2.08 


3000 


0.75 


5900 


0.78 


280 


3.5 


Lorentzian (7 = 0.5) 



M* = 2.5M , L* = 22L Q , d = 160 pc, A K = 2.0, i = 60°, R in = li?*, r = 5i? 

t Models L12, L14, and L15 use LTE level populations. Models U2 and V2 use non-LTE 

level populations. 



Table 2. 1548c27 Model Parameters 





Model 


R* 

(Rq) 


(K) 


Pc 


(K) 


Pi 


( g cm -2 


) 


<f>v 






Gl_10 
G5_9 


5.9 
5.9 


9800 
9800 


0.75 
0.75 


9700 
13000 


0.6 
0.7 


7100 
800 


2.5 
2.0 


Lorentzian 
Lorentzian 


(7 
(7 


= 1) 
= 5) 



M* = 4.OM , L* = 2OOL , d = 2400pc, A K = 0.62, i = 30°, R in = 1R*, r = 5.9i? 
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Figure Captions 

Figure 1. The v=2-0 bandhead emission from WL 16 which shows the characteristic shape 
of bandhead emission from a rotating disk: a blue wing, a shoulder, and an intensity peak 
redward of the bandhead. The spectral region covered by our observations includes the 
R39 - R62 lines. 

Figure 2. Schematic depiction of the convolution of the double-horned rotational broaden- 
ing function (heavy line) with the rest distribution of lines near the bandhead (light line). 
The rotational broadening function is the profile of an isolated line from an inclined Keple- 
rian disk that has a monotonically decreasing intensity distribution between an inner and 
outer radius. The rest spectrum of lines near the bandhead includes line blending due to 
local line broadening sources and assumes LTE level populations at a characteristic inner 
disk temperature. 

Figure 3. Model fits to the v=2-0 bandhead region of WL 16 (solid line) superposed on 
the observed spectrum (histogram): (a-c) LTE Models L12, L14, and L15; (d) Model L12 
including the H~ line layer continuum opacity (see text); (e,f) non-LTE Models U2 and 
V2. The difference between the upper and lower solid lines in each panel indicates the 
model continuum deficit or excess relative to the observed continuum level. 

Figure 4. Relative LTE optical depths (see text) for representative 5-3 (light solid line), 
4-2 (dotted line), 3-1 (short dashed line), 2-0 (long dashed line) lines in the 5-3 bandhead 
region compared with that for a line at the 2-0 bandhead (heavy solid line). 

Figure 5. The u=2-0 and u=5-3 bandhead regions of the spectrum of 1548C27. 

Figure 6. Non-LTE model fits to both the u=2-0 and u=5-3 bandhead regions of 1548C27 
with (a,b) Model Gl_10 (7 = 1 kms" 1 ) and (c,d) Model G5_9 (7 = 5 kms" 1 ). The u=2-0 
bandhead region has been scaled by 95% to match the continuum in both spectral regions. 

Figure 7. (a,b) The spectral fit obtained with the parameters of Model Gl_10 and the 
(artificial) extension of LTE level populations down to 600 K (r = 100i?*). (c,d) The 
spectral fit obtained with the same model parameters, non-LTE level populations, and 
an abrupt truncation of the emission beyond the radius at which T = 1500 K. As in 
Figure 6, the u=2-0 bandhead region has been scaled by 95% to match the continuum in 
both spectral regions. When compared with the fits shown in Figure 6, these figures show 
that non-LTE level populations introduce a sharp outer radius to the emission, an effect 
that is easily diagnosed by fitting these bandhead regions. 
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